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The weak lensing (WL) distortions of distant galaxy images are sensitive to neutrino masses by 
probing the suppression effect on clustering strengths of total matter in large-scale structure. We use 
the latest measurement of WL correlations, the CFHTLS data, to explore constraints on neutrino 
masses. We find that, while the WL data alone cannot place a stringent limit on neutrino masses due 
to parameter degeneracies, the constraint can be significantly improved when combined with other 
cosmological probes, the WMAP 5-year (WMAP5) data and the distance measurements of type-la 
supernovae (SNe) and baryon acoustic oscillations (BAO). The upper bounds on the sum of neu- 
trino masses are ^m^, = 1.1, 0.76 and 0.54 eV (95% CL) for WL-HWMAP5, WMAP5-hSNe-hBAO, 
and WL-fWMAP5-|-SNe-|-BAO, respectively, assuming a flat ACDM model with finite-mass neu- 
trinos. In deriving these constraints, our analysis includes the non-Gaussian covariances of the WL 
correlation functions to properly take into account significant correlations between different angles. 

PACS numbers: 04.50.-|-h, 98.80.Cq, 98.80.-k 

I. INTRODUCTION 

There are growing evidences that neutrinos have finite masses. The atmospheric, solar, reactor and accelerator 
neutrino oscillation experiments have confirmed that at least two of the three neutrino species have finite masses and 
their mass differences obtained from current experiments are lAm^^l ~ 7.7 x 10"^ eV^ and jAml^l ~ 2.4 x 10^"^ eV^ 
However, the absolute mass scale has yet to be determined because these experiments are only sensitive to the 
mass differences. Whereas kinematical probes such as tritium beta decay experiments and neutrinoless double beta 
decay experiments give upper bounds on the absolute neutrino mass 0, S 3 , cosmological observations provide a 
powerful, albeit indirect, means of constraining neutrino properties, indeed yielding a more stringent upper bound on 
the total neutrino mass at present. 

The finite-mass neutrinos affect structure formation via their effects on cosmic expansion history as well as on the 
evolution of perturbations in structure formation, which imprint characteristic signatures onto the cosmic microwave 
background (CMB) anisotropics and the mass clustering history in large-scale structure at low redshifts (see @ for a 
thorough review) . In particular the neutrinos slow down the growth of total matter clustering compared to a pure cold 
dark matter (CDM) model, because the neutrinos with large thermal velocities cannot be trapped by the gravitational 
potential well due to CDM inhomogeneities on small scales - the free-streaming effect Q. Combining the CMB 
information with low-redshift cosmological probes allows to efficiently track the suppression effect on gravitational 
clustering during the matter-radiation equality and low redshifts 0, Q- I^i fact, with the advent of high-precision 
cosmological probes, extensive efforts have been made in order to derive a stringent upper limit on neutrino masses 
by using/combining CMB and clustering properties of galaxies and Lyman-a forests d, 

[3, [m, E [il Q E H [13 ■ 

There is another vital probe of mass clustering: weak lensing (WL) or the so-called cosmic shear, the bending of light 
by intervening mass distribution that causes images of distant galaxies to be distorted (see [Tsl for a thorough review). 
These sheared source galaxies are mostly too weakly distorted to measure the effect on single galaxies, however, the 
lensing signals are measurable by correlating the different galaxy images - the shear correlation functions. The 
weak lensing probes the distribution of total matter along the line of sight, i.e. this method is free of galaxy bias 
uncertainties. Various groups have measured the weak len sing correlations and have also used the measurements to 
constrain cosmological models since the first detections [l9l . |20| [2ll . [2^ . Moreover there are several studies that assess 
the ability of future lensing surveys to constrain the neutrino masses [2^, [13, HE ■ However, the actual use of WL 
measurements for constraining the neutrino masses has not been explored yet. 
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In this paper, therefore, we pursue a new constraint on neutrino masses by using the latest WL data set, the 
Canada- France-Hawaii Telescope Legacy Survey (CFHTLS) [l^l, and also by combining it with other cosmological 
probes, the WMAP 5-year data, [1^, the type la supernova data [2§|, and the baryon acoustic oscillation experiment 
[30[. For this purpose, it is critically important to take into account the covariances of the WL correlation functions 
used, because there are significant cross-correlations between the WL correlations of different angles - the off-diagonal 
terms of the covariances are non-vanishing even for a pure Gaussian field. We use the halo model approach developed 
in (sTj to compute the covariances, and also discuss the effect of a possible uncertainty in the covariancc estimation 
on the final neutrino constraint. 

The paper is organized as follows. In Section II-A we define the WL correlation functions in terms of cosmological 
parameters and discuss how finite-mass neutrinos affect the WL signals. In Section II-B we describe our model to 
compute the covariances of the WL correlation functions. After giving a brief description of other cosmological probes 
(CMB, SNe and BAO) in Section II-D, we show the main results for the neutrino mass constraints in Section III. 
Section IV is devoted to conclusions. 



II. METHODOLOGY 



A. Preliminary 

In this section we briefly review how cosmic shear observable is related to cosmology, with particular attention to 
the effect of finite-mass neutrinos on cosmic shear. 

The cosmic shear field is the weighted mass distribution integrated along the line of sight. The cosmic shear 
fields are measurable only in a statistical sense. A most conventional method used in the literature is the two-point 
correlation based methods such as the two-point correlations of the shear field. The Fourier transformed counterpart 
is the shear power spectrum, P^, which is obtained by projecting the 3D matter power spectrum, Ps(k), weighted 
with the Icnsing kernel: 
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where x is the comoving distance along the light ray, xh is the distance to the Hubble horizon, i is the 2D wavevector 
perpendicular to the line of sight, i/o ~ lOO/i km s~^ Mpc""'^ is the Hubble constant, and fxix) is the comoving 
angular diameter distance out to a distance x- Note that in practice xh is taken to be the comoving distance out to 
a redshift where a source galaxy distribution is sufficiently decaying. The function n{z) is the redshift distribution of 
source galaxies, for which we adopt the following form characterized by three parameters according to (27| : 
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where r(a;) is the Gamma function. Note that n{z) is normalized so as to satisfy the condition J^dzn{z) = 1. The 
parameters a, (3, and z^ are calibrated by using the CFHT Deep Fields and the VIRMOS VLT Deep Survey [s^ . 
Because the uncertainties in parameters a and (3 are small they were fixed to 0.838 and 3.43, respectively, while the 
parameter Zg will be marginalized over assuming a Gaussian prior in the following analysis: Zg = 1.172 ifc 0.026, as was 
done in [3^. For I ^ 100 the major contribution to Pk{1) comes from nonlinear clustering (e.g., see Fig. 2 in [s^l). As 
described below, we employ the fitting formula in [s^ to compute the nonlinear Ps{k) from the input linear matter 
power spectrum. For the calculation of the linear matter spectrum, we use the CAMB code (36j up to the wavenumber 
piit — 25 Mpc ^, and small scale solutions of Hu & Eisenstein [s^l for k > fcgpiit- We chose this wavenumber 
because the bulk of information in WL comes from wavenumbers np to k/h ^ 1 Mpc~^ [1^. We confirmed that the 
results are stable against the change of /cgput. 

The formulae to map the linear matter power spectra to the nonlinear ones have been derived and tested against 
A'^-body simulations for CDM cosmologies. For a mixed dark matter model, the nonlinear matter power spectrum has 
yet to be understood, except for the initial attempts to study the weakly nonlinear regime of gravitational clustering 
based on the perturbation theory fs^ . |40| and the hybrid A'^-body simulation [ilj . In this paper we employ a similar 
prescription to that developed in |39l | (also see p6l |) in order to compute the nonlinear power spectrum for a mixed 
dark matter model. Our method is as follows. First, recall that, for a range of neutrino masses of interest, the 
finite mass neutrinos have large free-streaming scale below which the density perturbation of neutrinos is negligible; 
kfs ~ 0.04(mj//0.1 eV)(r2mo/0.27)^/^ hMpc~^ (see Appendix A in [1| for the details), where fimo is the present-day 




FIG. 1: Left panel: A fractional difference between the matter power spectra for a concordance ACDM model {Qmo ~ 0.3) 
and a model with finite-mass neutrinos {Qcho ~ 0.27 and fii,o — 0.03). Note that the total matter density r2mo(= ^cbo + ^vo) 
and other cosmological parameters are fixed for the two models. The solid curve shows the model prediction including the 
correction of nonlinear mass clustering (see the text for the details), while the dashed curve shows the linear-theory prediction. 
The finite-mass neutrinos cause a suppression in the power spectrum amplitudes on scales below the free-streaming scale. The 
suppression effect is enhanced over transition scales between linear and nonlinear regimes. Right panel: A similar plot, but for 
the lensing power spectrum as a function of multipoles I. 



where 



energy density of total matter in units of the critical density. Hence we assume that the neutrino density perturbations 
stay in the linear regime, which can be precisely computed by solving the hierarchical Boltzmann equations coupled 
to the Einstein equations using the CAMB. On the other hands, the nonlinear clustering is driven mainly by the 
density perturbations of CDM plus baryon. In this setting, the density perturbation of total matter (CDM, baryon 
and neutrinos), which generates the gravitational potential causing gravitational lensing via the Poisson equation, is 
given as 

Otot = - , _ — = /cbdcb + Ji^O^ , (3) 

Pc+ Ph+ 

_ ficO + f^bO r _ ^I'O c _ ^Pch + Spb , ^ 

Jcb = , J I- = t; — , Ocb = — — — • (.4) 

i 'mO S 'mO Pc + Ph 

The subscripts 'c', 'b', V and 'cb' stand for CDM, baryon, finite mass neutrinos, and CDM plus baryon, respectively, 
r^io (i = c, b and so on) denotes the present-day energy density of the i-th component, /, is its fractional contribution 
to the total energy density, and 5]^ represents the linear density perturbation of neutrinos. Note that /cb + /ly = 1- 
We neglect the gravitational dragging force on nonlinear density perturbations of baryon and CDM from neutrinos. 
As stated above, we use the mapping formula of [s^ to compute the nonlinear density perturbation of CDM plus 
baryon from the input linear one. Then the nonlinear power spectrum of total matter is given by 

= (O = fl {{5lf) + fl ((O^) + 2/./cb (SiSl^,) , (5) 

where we have introduced the superscript 'NL' to explicitly state the quantities in the nonlinear regime. The sensitivity 
of the nonlinear power spectrum depends to neutrino masses arises from the growth rate in each density perturbations 
as well as from the dependence of each terms in Eqn. (O on /j, . 

One may imagine a more crude approach that the nonlinear power spectrum is mapped by inserting the input linear 
power spectrum of total matter into the fitting formula treating the total matter as a single fluid. We have checked 
that the result obtained in this approach is not so largely different from our fiducial approach of Eqn. ^ as long as 
neutrino masses are small. Even so we believe that the method above is more sensible in a sense that the method 
includes the perturbation theory result [soj when the perturbations are in the weakly nonlinear regime. 

The effect of finite mass neutrinos on Pg'^^k) is depicted in the left panel of Fig. [2 showing the relative difference 
between the nonlinear power spectra with and without finite mass neutrinos, where f2mo is kept fixed. For comparison 
the corresponding result for the linear power spectra is also shown, where the denominator and numerator in APs/Ps 
are both the linear spectra. On very large distance scales (i.e. very small fc), well beyond the neutrino free-streaming 
scale, the neutrino effect is absent, as the neutrinos can cluster together with the CDM plus baryon perturbations. 



4 



For intermediate scales, corresponding to k ^ [0.01 — 0.1]hMpc~^ for the case with ~ 0.5 eV, the neutrinos cause a 
characteristic scale-dependent suppression in the matter power spectrum amplitude. Very interestingly, the neutrino 
suppression effect is enhanced on scales around k ~ l/tMpc~^ compared to the linear theory prediction, which is 
consistent with the perturbation theory result found in [39|. Perhaps even more surprisingly, the amount of the 
neutrino suppression effect becomes similar to that of the linear theory prediction on the smaller scales k ^ lhMpc~^ , 
well below the free-streaming scale. That is, the neutrino effect appears just as a constant offset in the overall amplitude 
on the small scales, even though the power spectrum amplitude itself is significantly boosted in the nonlinear regime 
compared with the linear theory. This seems consistent with the result indicated by the hybrid A^-body simulations 
[4l| . although the simulation resolution may be insufficient to follow this nonlinear regime. The amount of neutrino 
suppression on scales down to nonlinear regime is roughly given by APg/ Pg ~ — 8/^ which is sensitive to total neutrino 
mass as /j, ~ ^ ?7ii,/[nino/i^(94.1 eV)], but not sensitive to the mass scale of individual neutrinos. 

There is a good rationale to believe that our modeling of the nonlinear power spectrum is fairly reasonable according 
to the previous findings in the literature. As accounted for in the fitting formula [35|, the nonlinear spectrum amplitude 
at a given scale is very sensitive to the local amplitude as well as the local spectral slope rioff = c? In P}' /d\nk of the 
input linear power spectrum at that scale, where the power spectrum is a decreasing function with increasing k at 
relevant nonlinear scales (also see [12] )• In particular a change in the local spectral slope induces a scale-dependent 
modification in the shape of the nonlinear power spectrum over a range of the transition scales between the linear and 
nonlinear regimes (also see Fig. 1 in [4^). Having in mind these properties, the result in Fig. [1] for the input linear 
power spectra explicitly shows that including finite-mass neutrinos leads to a change in the local spectral slope over a 
range of the intermediate scales, causing the enhanced neutrino effect in the weakly nonlinear regime. On very small 
scales, the local spectral slope is unchanged, therefore, the neutrino effect becomes scale-independent. 

Thus the finite mass neutrinos affect the lensing power spectrum via the effect on the matter power spectrum, even 
for a fixed r^mOj which keeps the lensing efficiency function unchanged. The right panel of Fig. [1] shows the neutrino 
effect on the lensing spectra. The neutrinos suppress the lensing spectrum amplitudes over all the angular scales 
probed by the current data. Note that the baryon oscillation wiggles are completely washed out due to the projection. 
Thus, since the shear correlation function is a Fourier-transform of the lensing power spectrum, the neutrino effect can 
be potentially extracted from the measured shear correlations. In particular, if we can use the small angular scales 
with a robust model prediction, the accuracy of constraining neutrino masses can be significantly improved in the 
nonlinear regime compared to the linear regime, due to the increased signal-to-noisc ratio by the boosted amplitudes 
of shear correlations. 



B. The CFHTLS result for cosmic shear measurement 



The data we will use to put limits on neutrino masses is taken from the latest Canada-France-Hawaii Telescope 
Legacy Survey (CFHTLS or we will often call WL for simplicity) [l^. It covers over 57 square degrees, which drops 
to the effective area of 34.2 square degrees after taking into account the masking corrections of the surveyed region, 
and comprises about 1.7 x 10^ galaxies whose shapes were used for the lensing shear correlation analysis. In their 
paper, the shear correlation is given in terms of correlation function ^, 7, and Map, all of which are related to the same 
convergence power spectrum (Eqn. [T]) with different filtering. They found that these three statistics give consistent 
results on cosmological parameter estimations. In this paper we therefore use one of the correlation functions, ^. 

It is useful to decompose the cosmic shear correlations into the E (irrotational) and B (rotational) mode contribu- 
tions, because the cosmological shear field creates the i?-mode field alone. The i?-modc correlation is related to 
the convergence power spectrum P„ [il, IH, \^ as 

UO) = —P^(i)j„{i9), (6) 

where 6 is the separation angle between galaxy pairs, and Jg is the zero-th order Bessel function. In practice, since the 
shear power spectrum is difficult to measure due to the complex survey geometry, the E'-mode correlation function is 
estimated from the measured correlation functions ^± (0) as proposed in [43j : 



e+(^?)+e-(e)+ —i-{B')[A-l2-^ 



(7) 



where £,±{0) is the two-point correlations of ellipticities between paired galaxies separated by angle 6*; the compo- 
nent is the tangential or radial ellipticity component with respect to the line connecting the two galaxies, while the 
component 'x' is its 45 degree rotated component. Thus, the transformation relation of S^e from the measured ^±(6*) 
includes an integration of ^_ up to an infinite separation 6 —>■ 00. However, in practice, ^e needs to be estimated from 
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FIG. 2: The data points show the measured shear correlation functions, (,e{9), at each angular bins, which are taken from 
the CFHTLS result in [23]. The error bars around each data points are computed from diagonal terms in the inverse of the 
covariance matrix that includes contributions from the shot noise of intrinsic galaxy ellipticities and the Gaussian and non- 
Gaussian sample variances (see the text for the details). The solid curve is the model prediction for the ACDM model with 
finite- mass neutrinos which best matches the WL measurement. The dotted curve is the best-fitting model prediction for the 
joint fitting of WL-|-WMAP5-l-SNe-l-BAO as will be shown below. Note that the best-fitting model has the total neutrino mass 
of 'Y^rriv = 0.25 eV. To demonstrate the effect of finite-mass neutrinos on ^_b, the dashed curve shows the model prediction 
where the neutrino mass is changed to J]] = 0.54 eV, roughly at two sigma upper bound for the joint fitting, and other 
cosmological parameters are fixed to their best-fitting values. 

a finite range integration that may introduce a residual uncertainty in the estimated S^e{0). We ignore this uncertainty 
in our paper for simplicity. 

Fig. [2] shows the measured shear correlation function taken from Table B.l in [l^], comparing with the model 
predictions which best fit to to the CFHT WL data and the case combined with the WMAP5 and other distance 
measurements (see below for the details). The error bars at each angular bins are computed from the diagonal 
components of the inverse of the covariance matrix, and the data points between different separation angles are 
significantly correlated to each other as described below in detail. 

C. Covariances of shear correlation functions 

The covariance matrix of shear correlation function quantifies (1) statistical uncertainties in measuring the shear 
correlation function at each angular bin, and (2) how the correlation functions of different angular bins are correlated 
with each other. Following the method in [3l|, |43| (more explicitly see Eqn. [29] in [3l|), the covariance matrix of the 
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FIG. 3: The projection matrix \Sim\ for the principal component decomposition of the normahzed covariance matrix is pfotted, 
for first three eigenmodes that have largest differential contributions to the cumulative signal-to-noise ratio, {S/N), defined by 
Eqn.(l9ll. The left panel shows the results obtained when using the halo model developed in [3l|] to compute the covariance 
matrix. The total signal-to-noise ratio S/N — fO.6. The right panel shows the result when the fitting formula in Semboloni 
et al. (2007) [31 is employed, yielding (S/N) = 8.89. For the real-space correlation function, the projection matrix for each 
eigenmodes has a broad tail, refiecting significant correlations between different angles. It is also found that the contributions 
from each separation angles become different depending on which model of the covariance matrix is used. 



shear correlation function, C, can be expressed as 

[C],, EE COY[UO^),U0,)] 

= (5/;covsN + ^y^ idiuieMie.j)p^{i) + ^-^ idi i'di'Mi9,)Mio,)f,{ij'), (8) 

where Sfj denotes the Kroncckcr-typc delta function defined so that = 1 when 9i = 9j within the bin width, 
otherwise 0, As is the survey area in units of steradian and is the angle averaged trispectrum of convergence. Note 
that the survey area is set to Ag — 34.2 square degrees as described above, and we have ignored the boundary and 
geometry effects of surveyed region for simplicity. 

The first term on the r.h.s. of Eqn. ([8]), Sf^ CovsN, denotes the the shot noise contamination due to intrinsic galaxy 
ellipticities, contributing only to the diagonal terms of the covariance matrix with 9i = 9j . The shot noise is determined 
by the rms intrinsic ellipticities and the total number of galaxy pairs that are available from the surveyed region in 
the separation angle of a given bin width. The shot noise can be estimated directly from the data: for example, 
since the coherent shear signals can be erased by randomly rotating each galaxy images by arbitrary angles, the shot 
noise contamination can be estimated from variations in the correlation functions that are repeatedly measured after 
randomization of galaxy image orientations. Our analysis uses the shot noise at each bin given in the column labeled 
as "JCb" in Table Bl of [13. 

The second and third terms in Eqn. ([5]) denote the Gaussian and non-Gaussian sample variances, respectively, 
arising due to the imperfect sampling of shear correlations from a finite survey region. Since the lensing correlations 
probe the mass distribution in large-scale structure along the line of sight, modeling the sample variances requires a 
knowledge about the power spectrum and trispectrum of the 3D mass distribution, where non-Gaussian errors arise 
once the mass distribution probed by lensing resides in the nonlinear clustering regime. Thus the sample variances 
depend on the mass clustering strengths, i.e. on the underlying cosmology. There are several important features of 
these sample variances. First, contrary to the power spectrum covariance, the Gaussian term is non-vanishing for the 
off-diagonal components of the covariance when 9i ^ 9j : there are always correlations between the shear correlations 
of different angles. Second, the non-Gaussian sample variance contributes to both the diagonal and off-diagonal terms 
of the covariance matrix. The non-Gaussian errors become more significant on smaller angular scales where the lensing 
signals are more affected by the nonlinear regime. 

Thus taking into account the covariance is critically important in order not to have too optimistic parameter 
estimation from the measured correlation functions. However, only a few previous works have studied the covariances 
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of cosmic shear correlations, based on ray-tracing simulations [i^ and analytic methods [3l|, [13, [3 . In particular, 
the importance of non-Gaussian errors is not fully understood yet. In this paper, to estimate how an uncertainty in 
the covariance affects the results, wc employ two models of the non-Gaussian covariances: 

(A) Compute the covariance matrix using the dark matter halo approach developed in [3l[. 

(B) Compute the non-Gaussian error contribution to the covariance by multiplying the fitting formula derived in 
p6l | with each matrix elements of the Gaussian sample variance (the second term in Eqn.[8]). 

For the first approach we also include sources of new non-Gaussian errors that inevitably arise for a finite survey 
area H^, although the new contribution seems not so significant as studied in [3l|. Model A is our fiducial model to 
compute the covariance, but we will study how our results change when using Model B for the covariance evaluation. To 
compute the covariance we need to specify cosmological model, and our covariance evaluations assume the concordance 
ACDM model which is consistent with WMAP5. 

It is useful to study how the two-point correlation functions of different separation angles are correlated with each 
other. To do this, we study the cumulative signal-to-noise (S/N) ratio for measuring the shear correlation functions. 
Using the normalized covariance matrix defined as = Cij /[£,E{(^i)^E{()j)], the S/N is defined as 

^ ^ — l m—l \^— 1 / 

where C is the inverse matrix of C and Nb is the number of separation bins {Nb ~ 34 for CFHTLS as shown in 
Fig. [21). In the second equality we have used the principle component decomposition given as Cij = X^mLi Sim^mSjm, 
where Am is the m-th eigenmode and the projection matrix Sim describes how the i-th covariance element contributes 
to the TO-th eigenmode. Note that the inverse matrix of C is given by [C ]ij = J2m''=i ^im{^/^m)Sjm, and the 
matrix Sim satisfies the conditions S = S"^ ,J2m''=i ^imSjm = ^ij, and X^mLil'^'im.)'^ = 1- Eqn. ^ expresses the S/N 
as a sum of the contributions from independent eigenmodes. 

Fig-Elshows the projection matrix of first three eigenmodes that have largest differential contributions to the [S/NY 
in Eqn. ([9]), using cither of the two models above to compute the covariance matrix. It is worth noting that Model A 
gives S/N = 10.6, while Model B gives S/N = 8.89, because Model B predicts more significant non-Gaussian errors 
than Model A as studied in [sij. For both the models, each projection matrix gives contributions over a wide range 
of separation angles, reflecting signiflcant correlations between different angles. Thus there are fewer independent 
modes than the 34 bins of separation angles in Fig. [2| for both these two cases, first 15 eigenmodes give about 95% 
contribution to the (S/N)^. Also note that the angular scales mainly contributing to the total S/N are different in 
these two models: for Model B based on the method in Semboloni et al. (2007), the small angular scales are more 
important in constraining cosmology. 

We will use the log likelihood analysis to estimate cosmological parameter: 

where ^'^^ and ^j-e the measured correlation function and the model prediction, respectively, and is the 

inverse of the covariance matrix given as [C~^]ij = £,Eidi)[C ]ij£,E{Sj)- 



D. Likelihood analysis 

Our likelihood analysis includes a fairly broad range of cosmological parameters that affect a ACDM cosmology 
with finite mass neutrinos. We will work on 8 parameters given as 

p= (^7bo/^^^^dmo/^^^,eA,T,ln(10l°A2,),n„z,) , (11) 

where 9a is the so-called acoustic scale related to the distance to the last scattering surface [s^l, r^dmo^^ = {^co + 
^i^o)h^ stands for the dark matter density today, and and ri, are the amplitude and spectral index of the 
primordial curvature power spectrum defined at A: = 0.002 Mpc~^ |ll| . respectively. The parameter space we explore 
in the following is: fly,ah^ ^ [0.005,0.1], f^dmo^^ = [0.01,0.99], = [0,0.9], Oa = [0.3,10], r = [0.01, 0.8], In lO^^A^ = 
[0.5, 6.0], Ug — [0.5, 1.5], Zg — [0.9, 1.4]. Throughout this paper we assume the standard three flavor neutrinos that are 
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FIG. 4: The contours show the marginalized constraints (68% and 95% CL) for the m^, nnio)-subspace, obtained by fitting 
the WMAP5 data to the ACDM model with finite-mass neutrinos. 

completely degenerate in mass. Note that the parameter Zs specifies the galaxy redshift distribution that is needed 
to predict the shear correlation functions as described below Eqn. ([2]). and docs not affect other cosmological probes 
such as WMAP5, SNe and BAO. Assuming flat priors for the cosmological parameters except Zg, we employ the 
MCMC method [51| to explore cosmological parameter estimations in the multi-dimensional parameter space given 
cosmological observables. For each run four parallel chains were computed and the convergence test was made based 
on the Gelman and Rubin statistics ("i?— 1" statistics ^^). Our chains have typically 40000 points and i? — 1 ^ 0.05. 
In what follows we mention other cosmological probes used in our likelihood analysis, which adds complementary 
information to the cosmic shear data. 



1. CMB 

The CMB anisotropics are generated mainly until epochs before the last scattering surface of the decoupling epoch, 
z ~ 1089. Hence, if neutrino species are massive enough as m,j > 0.5 eV (or Y^m^ > 1.5 eV), they became non- 
relativistic before the recombination epoch In this case, the finite-mass neutrino causes a significant effect on 
the CMB spectrum mainly in two ways. First, for a fixed energy density of total matter (i.e. fJmO = ^cbo + ^i^o 
fixed), the matter-radiation equality is delayed by the presence of finite-mass neutrinos. Hence this makes the heights 
of acoustic peaks higher, especially that of the first peak, in the CMB spectrum due to a greater early integrated 
Sachs- Wolfe effect. Second, the neutrinos cause an additional decay of the gravitational potential on small scales. This 
effect boosts the amplitudes of acoustic oscillations at higher multipoles through the dilation effect due to diminishing 
gravitational potential [53|. Through these effects, a precise measurement of the CMB spectrum can be used to place 
a limit on neutrino masses. In this paper we use the WMAP 5-year (WMAP5) data, and explore the likelihood which 
can be evaluated by the code publicly available at LAMBDA web site (63 |. 

For the masses below rrii, < 0.5 eV (X)"*!^ S l-5eV), the neutrinos affect the CMB spectrum mainly through the 
effect on the angular diameter distance out to the last-scattering surface. In this case the effect is degenerate with 
other cosmological parameters such as fimo and h 0]. Thus other cosmological probes complementary to CMB are 
needed to break the parameter degeneracies in order to explore the small mass scales of neutrinos. 
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2. Type la SNe 

The distant type la supernovae (SNe) are one of such probes. The use of SNe as a standard candle aUows us to 
estimate the luminosity distance - redshift relation. The estimated distance, more precisely the distance modulus, 
can be compared with a model prediction for a flat ACDM model with finite-mass neutrino contribution: 

DLiz) = '^^^-r^ r dz' [ao(l + zT + f^cbo(l + z'f + n^z') + Qao] , (12) 

where n^o and f2Ao sre the energy densities of photon and cosmological constant, respectively, f2Ao ~ 1 ~ i^mo ^^i a 
flat universe, and n,y{z) is the energy density at a redshift z. The contribution from radiation is negligible for the 
luminosity distance out to SNe at z < 2. In addition, since the current cosmological probes are only sensitive to 
neutrino masses greater than ^ O.leV as will be shown below, neutrinos in this mass scale behave as non-relativistic 
matter over the redshift range of SNe data. Therefore, if we restrict our attention to a flat ACDM universe, the 
measured distance modulus of SNe constrain the information mainly on the total matter density Qrno- We use the 
latest compilation of SNe data set given by [2^. The SNe catalog provided consists of 307 SNe and is carefully 
calibrated taking into account various sources of the systematic errors. 



3. BAO 



Baryon acoustic oscillation (BAO) offers another probe. CMB measures the acoustic oscillations in the primordial 
photon-baryon plasma, which can be used to derive the angular diameter distance to the last scattering epoch, 
z ~ 1089. The same oscillation feature is imprinted onto the late-time matter power spectrum probed by a massive 
galaxy survey [H, [11]. This baryon acoustic oscillation feature therefore provides a standard ruler by which one 
can measure the distance to the redshift where the bulk of galaxies are observed. Eisenstein et al. [301 analyzed 
galaxies from Sloan Digital Sky Survey (SDSS) and provided their result in a form of an effective distance measure 
by averaging the distances in the radial and transverse directions: 



Dv{z) 



{l + zfDl{z) 



-,1/3 



H{z) 



(13) 



where Da{z) is the angular diameter distance to the redshift z = 0.35. In this paper we use the following constraint 
on the distance parameter A provided in [soj : 



A[z = 0.35) EE £1^(0.35) 



0.35c 



0.469 



f— ) 
Vo.98/ 



-0.35 



±0.017 



(14) 



Notice that the value of A is primarily sensitive to fJmo for a flat universe case. 



III. RESULT AND DISCUSSION 



A. Parameter constraints 



The WMAP team reported an upper limit on the total neutrino mass, X] "^j^ < 1.3 eV (95%CL), for a flat ACDM 
model plus one additional parameter of the total neutrino mass (TT| . Using the same likelihood function given by 
WMAP team, we find a consistent but slightly tighter limit, < 1.2 eV. This difference would be attributed to 

differences in the treatment of cosmological parameters as well as the priors between ours and WMAP5. In fact, we 
have found the consistent result if the same parametrization is adopted as theirs. In any case, these upper bounds 
are fairly close to the critical value ~ 1.5 eV as discussed in the previous section. 

Fig. [4] shows the results we obtained from the WMAP5 data: the inner and outer contours show the 68% and 
95% confidence level regions in the (ri,no, X)™i^)"Siibspace, marginalized over other parameters. The contours show 
a strong degeneracy between these two parameters, reflecting that the CMB constraints come mainly from the effect 
of finite-mass neutrinos on the matter-radiation equality. The delay in the matter-radiation equality caused by the 
presence of massive neutrinos is compensated by adding more non-relativistic matter at present, i.e. increasing r^mo- 

The results in Fig. [3] are not so encouraging from the perspective of complementarity between CMB and WL, 
because the WL constraints on the parameters fimo and ^ are expected to show a similar degeneracy curve. 
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FIG. 5: The improvement in parameter constraints for (fimo, X] ) obtained by combining the CFHT WL data with WMAP5 
(green-color contours), in comparison with the constraints for WMAP5 alone (black). The accuracy of flmo determination is 
improved by adding the WL constraint, because the WL amplitude is sensitive to SlmO as can be found from Eqn. ([T}. However, 
the constraint on Sm,y remains almost unchanged due to parameter degeneracies in the WL information. 

The finite-mass neutrinos cause suppression in the amplitudes of cosmic shear correlations on relevant angular scales, 
however, this effect is compensated by increasing fimo and/or the power spectrum normalization (see Fig. [T] and 
Eqn. [1]). Even so, since the degeneracy curves of CMB and WL constraints are not elongated to exactly the same 
direction, combining these two can always improve the constraints to some extent if they are consistent with each 
other. 

Fig. shows the constraints on the (Omo, X] ™i')"Pl^^s obtained when combining WMAP5 with the CFHT WL 
constraints. As expected from the above argument, the constraint is only slightly improved when the WL data is 
included. The upper bound on the total neutrino mass is ^rrii, ^ 1.1 eV for WL-fWMAP5 compared to '^i^ ^ 1-2 
eV for WMAP5 alone. The reason why the constraint is not significantly improved is also attributed to the quality of 
WL data: the current WL data with ~ 30 square degrees is not so powerful yet compared to the WMAP5 data. Finally, 
the weak lensing data itself suffers from a strong degeneracy between J7mo a-nd the power spectrum normalization, 
e.g. the flino-cs degeneracy, as shown in [27| : adding the WL data cannot yet much break the degeneracy between 
and flmo in the WMAP5 constraints. 

Since the geometrical probes such as SNe and BAO are sensitive to ilnM for a fiat ACDM model, by adding these 
probes into the WMAP5 and WL data, the constraint on neutrino masses can be further improved. In fact the 
WMAP team nicely demonstrated that combining SNe and BAO with the WMAP5 information tightens an upper 
limit on the neutrino masses by a factor of 2, resulting in ^mi, < 0.67 eV [T]| . 

Our analysis also confirmed this trend, but gives a slightly less improved limit, ^ < 0.76 eV, as shown in Fig. [B] 
The difference would be due to the fact that we are using different BAO constraint from that used by the WMAP 
team. The BAO constraint employed here was one derived in Eisenstein ct al. [soj using the spectroscopic sample 
of the SDSS Luminous Red Galaxies (LRGs) that probes the distance out to z = 0.35. This BAO measurement has 
a less constraining power than that more recently derived by Percival et al. [s^, which was used in the WMAP5 
analysis. Since they used a bigger galaxy sample combining the SDSS LRGs with the SDSS main galaxy sample and 
the 2dF galaxy sample, which probe the distances to redshifts z = 0.2 and z = 0.35. The constraints from Percival et 
al. are more powerful than those from Eisenstein et al.; for example, Percival et al. constrain fimo with the precision 
of Omo = 0.249 ±0.018 for a ACDM model, while Eisenstein et al. give = 0.273 ±0.025. In addition, as indicated 
from Fig. [5l the smaller best-fitting value of finio prefers a smaller J2 ''^f when the WMAP5 and BAO is combined. 
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However, a marginal tension between the two data sets of BAO has been discussed in the hteraturc, which leads to a 
possible inconsistency for a ACDM model (e.g., see Fig. 12 in [HI). This tension is not found if the BAO distance of 
z = 0.2 is not included. For these reasons, we employ the BAO constraint by Eisenstein et al. to derive a conservative 
constraint on neutrino masses. 

Fig. [6] shows how the neutrino mass constraints are improved when some or all of WMAP5, WL, BAO and SNe 
constraints are combined. Interestingly, adding the WL data into WMAP5+BA0+SNe does improve the constraint 
on the total neutrino mass, leading to the upper bound "^rrii, < 0.54 eV, because the degeneracy with flmo is more 
efficiently broken. It is also worth noting that a most significant degeneracy inherent in the weak lensing constraints, 
the riino~cr8 degeneracy, is efficiently broken by combining all the probes, as shown in Fig. [71 Our final results are 
^mo = 0.294 ± 0.02 and erg = 0.758 ± 0.03, improved by a factor of 10 and 5 from the constraints by the WL data 
alone. 

The shear correlation function for our best-fitting model is shown by the dashed curve in Fig. [21 in comparison with 
the measurement of WL and the model prediction (solid curve) that matches the WL data alone. The prediction 
of the best-fitting model to all the data combined (dotted curve) has systematically smaller amplitudes than the 
measurement data points, especially on large angular scales. However, it is worth noting that the change in 
between the two best-fitting models for the joint fitting and the WL data alone is only A^vvl — 2 due to significant 
correlations between the different data bins as implied in Fig. [3l Adding more massive neutrinos to the model with 
other parameters being fixed leads to further smaller WL amplitudes compared to the data. The dashed curve shows 
the result for such a model assuming ^ = 0.54 eV, a mass scale at about 2a upper limit for the joint fitting, 
whose value differs from the dotted curve by ^^Xwh — 2- 

The marginalized constraints and likelihood distribution of individual parameters of interest are summarized in 
Table Hand Fig. [SI respectively The precisions of constraining some parameters such as rimOj f^dmo^^j and ^ nii, 
are significantly improved when all the probes arc combined. 



TABLE I: Summary of the constraints on cosmological parameters and the marginalized errors for vACDM model 



Parameter 


WL 


WL±WMAP5 


WMAP5±SN±BA0 


ALL 




0.22 ±0.09 


0.12 ±0.0065 


0.11 ±0.0038 


0.11 ±0.0030 


Ho [km s"^ Mpc-i] 


74± 15 


65 ± 3.9 


68 ± 2.0 


68 ±1.9 




1.0 ±0.26 


0.95 ± 0.016 


0.96 ±0.014 


0.96 ±0.014 


ln(10"^A2,) 


2.8 ± 1.0 


3.2 ±0.046 


3.2 ± 0.040 


3.2 ±0.039 


Em. [eV] 


< 8.1 (95% CL) 


< 1.1 


< 0.76 


< 0.54 




0.64 ±0.15 


0.73 ± 0.045 


0.73 ± 0.059 


0.76 ± 0.033 




0.52 ±0.20 


0.33 ± 0.056 


0.29 ±0.021 


0.29 ±0.020 



B. Discussion of systematic errors 

There are several sources of systematic errors to affect weak lensing measurements. One of those is uncertainties 
in the source galaxy redshifts. As described in detail in [23|, the redshift distribution of source galaxies is carefully 
calibrated by using secure photometric redshifts of the CFHT deep fields, in combination with the VIMORS spectro- 
scopic sample, assuming that the deep fields contain a representative galaxy sarnple of the WL galaxies. We below 
estimate how a possible residual error in the source galaxy redshifts, reported in |27j |. affects our neutrino constraint. 
Since the amplitude of cosmic shear correlations is sensitive to mean redshift of galaxies, a most relevant parameter in 
the source galaxy distribution is the parameter in Eqn. ^ for which we have so far assumed Zg = 1.172 ±0.026 for 
the fiducial value and the la uncertainty. Note that we have also employed the Gaussian prior a{zs) = 0.026 around 
the central value in our MCMC analysis. 

The left panel of Fig. [Qlshows how the constraints on J2 ''^v ^'^^ '^8 change for the WL data combined with WMAP5 
if we put a more restrictive prior on Zs by a factor of 2 or 4 than the fiducial prior. The effects appear to be small: 
the confidence regions shrink only slightly. This is simply because the uncertainty in Za considered here is already 
sufficiently small compared with the statistical uncertainties in the shear correlation measurement. 

The right panel of Fig. [9l shows how a shift in the central value of Zg causes a bias in the best- fitting cosmological 
parameters, assuming the fiducial prior a{zs) = 0.026. The smaller or larger source redshift gives larger or smaller 
best- fitting values of ag, respectively (e.g., see [l^, [131), while the neutrino constraint is little affected. The effects 
of uncertainties in these parameters on cosmological constraints are found to be very small (~ 1%) in the present 
analysis. However, note that it will be of critical importance to have a precise knowledge on source galaxy redshifts 
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FIG. 6: Two-dimensional marginalized constraints on the neutrino mass and ^ImO are significantly improved when the 

constraints for WMAP5-I-WL are combined with the information from the geometrical probes, SNe and BAO. The black, 
green and red contom-s show the resuhs for WMAP5, WMAPS+SNe+BAO, and WL+WMAFS+SNe+BAO, respectively. The 
marginalized upper bound on the neutrino mass is given as ^ 0.54 eV for all the probes combined. 



especially when we have a sufficiently accurate measurement of cosmic shear that is more powerful in constraining 
neutrino masses than the CMB information [26j . 

As described in § III C[ there is an uncertainty in estimating the sampling variance contributions to the covariance 
matrix of the shear correlation functions. We have so far used the halo model approach developed in [3l| to compute 
the covariance matrix. Fig. [TO] shows the constraints on ^ and flmo obtained when using the different model for 
the covariance matrix evaluation, Model B based on [i^. These models are different in computing the non-Gaussian 
error contributions to the covariance. The confidence regions for Model B arc enlarged compared to those for Model A. 
The best-fitting parameters are only slightly changed: '}2,m,j = 0.25 and 0.21 for Model A and B, respectively. Model 
B predicts more significant non-Gaussian errors on angular scales of 6* ^ 10 arcminutes than Model A, i.e. stronger 
correlations between the shear correlation functions of different angular bins. The correlations reduce independent 
modes of the angular scales to be useful for constraining cosmology. This explains larger confidence regions for Model 
B. Also, as shown in Fig. [31 different angular scales of the shear correlation functions contribute to the cosmological 
constraints in a different way between Model A and B, causing a slight difference between the best-fitting parameters. 
Even so, it can be found that, thanks to a wide coverage of angular scales for the CFHT, the results for the two 
covariance models are consistent, implying that the large angular scales where the Gaussian errors are dominant in 
the covariance mainly contribute to the constraints on cosmological parameters. 

There are other sources of systematic errors such as those inherent in the galaxy shape measurements. Estimating 
this effect is beyond the scope of our paper, so we do not discuss those effects here. 
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FIG. 7: Two-dimensional marginalized constraints on erg and flmo, for WL+WMAPS (green) and WL+WMAP5+SN+BA0 
(red), respectively, in comparison with the constraints from WL alone (black). The yellow contours show the region from 
WMAP5+SN+BA0 (without WL). 



IV. CONCLUSION AND DISCUSSION 



In this paper we have investigated how the latest cosmic shear data set, the CFHT weak lensing survey, can constrain 
neutrino masses. The weak lensing correlations can probe clustering strengths of total matter free of galaxy bias 
uncertainty, therefore, it allows us to explore the suppression signal in the weak lensing amplitudes caused by the free- 
streaming effect of finite-mass neutrinos. Since the measured shear correlation functions do not have strong features 
in the angular dependence, the shear amplitudes carry most of cosmological information. However, due to the limited 
information for the current data sets, the shear amplitudes integrated over a range of source galaxy redshifts are only 
available: the redshift dependence of the shear signals is not. For these reasons, constraining neutrino masses from the 
weak lensing data alone suffers from severe degeneracies among cosmological parameters such as ^ m^, Qn-iO and the 
power spectrum normalization all of which are sensitive to overall amplitudes of the shear signals. Yet, interestingly 
enough, the neutrino mass constraint can be significantly improved when the weak lensing constraint is combined with 
other cosmological probes such as WMAP5, SNe and BAO, because the combination of these data sets helps efficiently 
break parameter degeneracies. The upper limits on the total neutrino mass found in this work are summarized as 
Y^m^ < 1.1 eV (WL+WMAP5), 0.76 cV (WMAP5-hSNe-hBA0), 0.54 eV (WL-f-WMAP5-hSNe-hBA0) at 95% 
confidence level, respectively. Note that we could not find any strong evidence of the lower limit for neutrino masses, 
i.e. a detection of the neutrino mass. This might be attributed to insufficient statistical precisions of the current data 
set 30 square degrees). 

Throughout this paper we have assumed a minimal cosmological model: a concordance ACDM model plus one 
additional parameter of the total neutrino mass. Thus our approach derives the best-available constraint on the 
neutrino mass from the weak lensing data combined with other cosmological probes. The constraint will be relaxed 
if additional parameters arc included. In particular, dark energy equation of state parameter (w = Pdo/pdo) may 
be most intriguing to include, because a change in the equation of state from the cosmological constant {w = — 1) 
does change the cosmic shear correlations in the shape and amplitude 3j|. The cosmological constraints from the 
geometrical probes, SNe and BAO, are also relaxed if including w as a free parameter. Therefore the accuracy of the 
neutrino mass constraint is degraded if w is included as a free parameter. This is beyond the scope of this paper, and 
will be presented elsewhere. 

Various weak lensing surveys are being planned and proposed so as to have a much increased survey area by a 
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FIG. 8: One-dimensional marginalized likelihood on individual cosmological parameter for WMAP5 alone (solid curve), 
WL+WMAPS (blue dashed) and WL+WMAFS+BAO+SNe (red dotted), respectively. Some of parameters, flmOy^dmoh^ , 
as and "^m^, are significantly improved by combining all the probes, where those parameters are sensitive to the overall 
amplitudes of cosmic shear correlations. 



factor of 100 compared to the CFHT survey, corresponding to a survey with area more than a few thousands square 
degrees. These surveys will promise to drastically improve the statistical accuracies of the cosmic shear measurements 
by a factor of 10 at each angular bins. Our results of Table |T] imply that such a survey may allow a precision 
of constraining the neutrino masses at level a(^m,y) ~ 0.8 eV with the weak lensing data alone, if systematic 
errors inherent in the measurements are well under control. Combining with other probes allows to further improve 
the neutrino mass constraint as demonstrated in this paper. In addition, if multi-color information is available for 
future surveys as proposed, photometric redshift information is available to estimate the redshift distribution of source 
galaxies in a more reliable fashion. This additional redshift information is extremely useful in that it enables to recover 
redshift dependence of the lensing signals, the so-called lensing tomography, allowing to efficiently break parameter 
degeneracies such as ag-^mo degeneracy with the weak lensing data alone [34| . Therefore, if lensing tomography is 
available, weak lensing (combined with CMB and others) would potentially allow a detection of the neutrino masses. 

In addition to the weak lensing technique, several cosmological probes have been proposed to constrain neutrino 
masses at accuracies of 0(0. 1) eV, such as CMB lensing ^5^, 21cm experiments [H^, galaxy surveys Q, cluster number 
count (g^I, and so on. All of these techniques are essentially based on the same effects of neutrinos, i.e., neutrinos 
suppress the power of density fluctuations on small scales. Therefore, these experiments including weak lensing allow 
us to make cross-checks of the cosmological constraint and hopefully to detect the finite mass of neutrinos, resolving 
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FIG. 9: Effects of possible residual uncertainties in the mean redshift of weak lensing galaxies on the cosmological constraints 
on C^m^, as) for WL+WMAP5. The mean redshift is parametrized by Zs in Eqn. ((5)), and the parameter is calibrated as 
Zs = 1.172 ± 0.026 according to [s^]. We have so far assumed the Gaussian prior of a{zs) — 0.026 around the central value 
in the parameter estimations. Left panel: Shown is how the accuracy of the parameter estimation is changed if we assume a 
more restrict prior by a factor of 2 and 4 than the fiducial prior. The difference is indistinguishable, implying that this error is 
negligible compared to the measurement errors in shear correlations. Right panel: The effect of a shift in the central value of 
Zs on the parameter constraints: the green and red contours show the results obtained when the central value of Zs is changed 
by about la statistical uncertainty to Zs — 1.146 and 1.198, respectively. The difference is very small. 



the hierarchy of neutrino masses with cosmological data sets. 

When this paper was under completion, similar works for constraining the neutrino mass with the CFHT data 
combined with other cosmological probes were put forward by [6ll . [6^ . In [6l| . they found an upper bound on 
neutrino masses which is consistent with ours. However, they also reported a lower bound ^ > 0.03 eV which 
we could not find here. The difference may be attributed to different data set used for BAO and SNIa, and different 



covariance matrices from ours. In |62| they derived a constraint on neutrino masses from CFHT data without 
combining WMAP5 data. Therefore their obtained upper bound is slightly weaker than ours but it is a consistent 
result. After we submitted our paper, a similar paper [63| was posted on the arXiv, in which further data set such as 
SDSS LRG and the other CMB data were used to obtain the upper bound on neutrino masses. 
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